World Nematic Crystal Model of Gravity 
Explaining the Absence of Torsion 
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We attribute the gravitational interaction between sources of curvature to the world being a 
crystal which has undergone a quantum phase transition to a nematic phase by a condensation of 
dislocations. The model explains why spacetime has no observable torsion and predicts the existence 
of curvature sources in the form of world sheets, albeit with different high-energy properties than 
those of string models. 



I. INTRODUCTION 

Present-day string models of elementary particles are 
based on the assumption that relativistic physics will 
prevail at all energy scales and, moreover, show recur- 
rent particle spectra at arbitrary multiples of the Planck 
mass. Disappointed by the failure of these models 1] to 
explain correctly even the the low-lying excitations, and 
the apparent impossibility of ever observing the char- 
acteristic recurrences, an increasing number of theoreti- 
cians is beginning to suspect that God may have chosen 
a completely different extension of present-day Lorentz- 
invariant physics to extremely high energies This 
philosophy has been advocated by one of the authors 
(HK) for almost two decades. In 1987, he proposed a 
simple three-dimensional euclidean world crystal model 
of gravitation in which dislocations and disclinations rep- 
resent curvature and torsion in the geometry of spacetime 
H. A full theory of gravity with torsion based on this 
picture is pubhshed in the textbook (6| (see also 

The simple 1987 model had the somewhat unaesthetic 
feature that the crystal possessed only second-gradient 
elasticity to deliver the correct forces between the sources 
of curvature, which for an ordinary first-gradient elastic- 
ity grow linearly with the distance R and are thus con- 
fining. In this note we would like to point out that the 
correct l/i?-behavior can also be obtained in an ordinary 
first-gradient world crystal with first-gradient elasticity 
by assuming that the dislocations have proliferated. This 
explains also why the theory of general relativity requires 
only curvature for a correct description of gravitational 
forces, but no torsion. Such a state of the world crystal 
bears a close relationship with the nematic quantum liq- 
uid crystals of condensed matter physics, first suggested 
by Kivelson et al. 0], and believed to be of relevance 
both for the quantum Hall effect ^3 ^"^^ bigh-Tc su- 
perconductors Tl'l. 

Our model will be formulated as before in three eu- 
clidean dimensions, for simplicity. The generalization to 
four dimensions is straightforward. The elastic energy is 
expressed in terms of a material displacement field Ui (x) 



as 



where 



"yW = -[a,Uj(x) +9jMj(x)] 



(1.1) 



(1.2) 



is the strain tensor and /i, v are the elastic the shear 
moduli. The elastic energy goes to zero for infinite wave 
length since in this limit Ui(x) reduces to a pure trans- 
lation under which the energy of the system is invariant. 
The crystallization process causes a spontaneous break- 
down of the translational symmetry of the system. The 
elastic distortions describe the Nambu-Goldstone modes 
resulting from this symmetry breakdown. Note that so 
far the crystal has an extra longitudinal sound wave with 
a different velocity than the shear waves. 

A crystalline material always contains defects. In their 
presence, the elastic energy is 
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(1.3) 



where ufj is the so-called plastic strain tensor describing 
the defects. It is composed of an ensemble of lines with 
a dislocation density 

ckji = e,jkdjdkUi{^) = (5i(x; L){hi + eiqr^qXr). (1.4) 

and a disclination density 

Oil = ^ijkdj(j>li = (5i(x;i)r2/, (1.5) 

where bi and ili are the so-called Burgers and Franck vec- 
tors of the defects. The densities satisfy the conservation 
laws 



diaik = -et 



mn^mn ; 



d^i = 0. 



(1.6) 



Dislocation lines are either closed or they end in discli- 
nation lines, and disclination lines are closed. These are 
Bianchi identities of the defect system. 
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An important geometric quantity characterizing dis- 
location and disclination lines is the incompatibility or 
defect density 



(1.7) 



It can be decomposed into disclination and dislocation 
density as follows 0: 

r]ij{x.)^9^j{x.) + ^dra [emmajn(x) + (i^j)-ey„a„„(x)] . 

(1-8) 

This tensor is symmetric and conserved 

a,77y(x)-o, (1.9) 

again a Bianchi identity of the defect system. 

It is useful to separate from the dislocation density 
(|1.4() the contribution from the disclinations which causes 
the nonzero right-hand side of ()1.6|> . Thus we define a 
pure dislocation density 



a,j{x) = ay(x) - ay(x) 



which satisfies dia^j 



0. Accordingly, we split 



'7uW =<-(x)-f ry,^-(x), 



(1.10) 



(1.11) 



where 



^r, W=2 [em^«a5„(x) + (i^j) -ey„a^„(x)], (1.12) 

and the pure disclination part of the defect tensor looks 
like (|1.8|l . but with superscripts il on ry^ and a^j. 

The tensors , 9ij , and rjij are linearized versions 
of important geometric tensors in the Riemann-Cartan 
space of defects, a noneuclidean space with curvature and 
torsion. Such a space can be generated from a flat space 
by a plastic distortion, which is mathematically repre- 
sented by a nonholonomic mapping 0)111 ^ Xi+Ui{x). 
Such a mapping is nonintegrable. The displacement fields 
and their first derivatives fail to satisfy the Schwarz in- 
tegrability criterion: 



Part IV in Vol. II of where the full one-to-one corre- 
spondence between defect systems and Riemann-Cartan 
geometry is developed as well as a gravitational theory 
based on this analogy. 

Let us now show how linearized gravity emerges from 
the energy H1.3|l . For this we eliminate the jumping sur- 
faces in the defect gauge fields from the partition function 
by introducing conjugate variables and associated stress 
gauge fields. This is done by rewriting the elastic action 
of defect lines as 



d^3 



1 



, (1.14) 



where ly = A/2(A -I- /i) is Poisson's ratio, and forming 
the partition function, integrating the Boltzmann factor 
^-E/ksT Qygp (j^j^Ui, and summing over all jumping sur- 
faces S in the plastic fields. The integrals over Ui yield 
the conservation law dia-ij — 0. This can be enforced as 
a Bianchi identity by introducing a stress gauge field hij 
and writing ct^ = = eikiejmndkdmhn- The double 
curl on the right-hand side is recognized as the Einstein 
tensor in the geometric description of stresses, expressed 



in terms of a small deviation n hj 



-Si-j of the metric 



H] — gij ^Hj ___ 
from the flat-space form. Inserting Gy into (|1.14|) and 
using (|1.7|) . we can replace the energy in the partition 
function hy E = E^^"'^'' + E'^''^ where 



-E 



dof - 



in ( 
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(1.15) 



where the defect tensor H1.8|l has the decomposition 

(1.16) 
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The defects have also core energies which has been ig- 
nored so far. Adding these for the dislocations and ig- 
noring, for a moment, the disclination part of the defect 
density in l|1.16|l . we obtain 



Tidisl _ 



2 ^" 



(1.17) 



{d,d,~d,d,)u{x.) ^ 0, id,dj~d,d,)dkui{x) ^ 0.(1.13) 

The metric and the afhne connection of the geometry in 
the plastically distorted space are gij = 6ij + diUj + djUi 
and Tiji = didjui, respectively. The noncommutativity 
of the derivatives in front of ui (x) implies a nonzero tor- 
sion, the torsion tensor being Sijk = (Tijk — ^jik)/'2- The 
dislocation density is equal to = eikiSkij, 

The noncommutativity of the derivatives in front of 
9fcU;(x) implies a nonzero curvature. The disclination 
density 6ij is the Einstein tensor dij — Rji — ^gjiR of 
this Einstein-Cartan defect geometry. The tensor 77^, 
finally, is the Belinfante symmetric energy momentum 
tensor, which is defined in terms of the canonical energy- 
momentum tensor and the spin density by a relation just 
like H1.8|l . For more details on the geometric aspects see 



We now assume that the world crystal has undergone a 
transition to a condensed phase in which dislocations are 
condensed. To reach such a state, whose existence was 
conjectured for two-dimensional crystals in Ref. f^e 
model requires a modification by an additional rotational 
energy, as shown in ^13] and verified by Monte Carlo sim- 
ulations in The three-dimensional extension of the 
model is described in . 

The condensed phase is described by a partition func- 
tion in which the discrete sum over the pure dislocation 
densities in is approximated by an ordinary func- 
tional integral. This has been shown in Ref. Q- The 
general integration rule is 



d^l6id-l)e 



1/20 



(1.18) 
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where has the components an = —ieijkdjUk/ \J —d^ . 
The Boltzmann factor resulting in this way from E'^'^^'^^^ 
plus (|1.17l) has now the energy 
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(1-19) 

The second term implies a Meissner-like screening of the 
initially confining gravitational forces between the discli- 
nation part of the defect tensor to Newton- like forces. 
For distances longer than the Planck scale, we may ig- 
nore the stress term and find the effective gravitational 
action for the disclination part of the defect tensor: 
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(1.20) 



A path integral over hij and a sum over all line ensem- 
bles applied to the Boltzmann factor e~^^^ is a simple 
Euclidean model of pure quantum gravity. The line fluc- 
tuations of ry^- describe a fluctuating Riemann geome- 
try perforated by a grand-canonical ensemble arbitrar- 
ily shaped lines of curvature. As long as the loops are 
small they merely renormalize the first term in the en- 
ergy (|1.20|l . Such effects were calculates in closely related 
theories in great detail in Ref. . They also give rise to 
post-Newtonian terms in the above linearized description 
of the Riemann space. 

We may now add matter to this gravitational environ- 
ment. It is coupled by the usual Einstein interaction 



E'' 



(1.21) 



where T'-' is the symmetric Belinfante energy momentum 
tensor of matter. Inserting for Gij the double-curl of 
hij we see that the energy (ll.2Uf) produces the correct 
Newton law if the core energy is Ec — SttG, where G is 
Newton's constant. 

Note that the condensation process of dislocations has 
led to a pure Riemann space without torsion. Just as 
a molten crystal shows residues of the original crystal 
structure only at molecular distances, remnants of the 
initial torsion could be observed only near the Planck 
scale. This explains why present-day general relativity 
requires only a Riemann space, not a Riemann-Cartan 
space. 

In the non-relativistic context, a dislocation conden- 
sate is characteristic for a nematic liquid crystal, whose 
order is translationally invariant, but breaks rotational 
symmetry (see 0, in the two dimensions and in 
the 2-|-l-dimensional quantum theory). The Burgers vec- 
tor of a dislocation is a vectorial topological charge, and 
nematic order may be viewed as an ordering of the Burg- 
ers vectors in the dislocation condensate. Such a mani- 
fest nematic order would break the low energy Lorentz- 
invariance of space-time. We may, however, imagine that 
the stiffness of the directional field of Burgers vectors is 



so low that, by the criterion of Ref. U6], they have under- 
gone a Heisenberg-type of phase transition into a direc- 
tionally disordered phase in an environment with only a 
few disclinations. In three dimensions, dislocations (and 
disclinations) are line-like. This has the pleasant conse- 
quence, that they can be described by the disorder field 
theories developed in ^3 in which the proliferation of 
disclinations follows the typical Ginzburg-Landau pat- 
tern of the field expectation acquiring a nonzero expec- 
tation value. A cubic interaction becomes isotropic in 
the continuum limit [l^ (this is the famous fiuctuation- 
induced symmetry restoration of the Heisenberg fixed 
point in a (p^-iheory with 0(3)-symmetric plus cubic in- 



teractions jl9j . The isotropic phase is similar to what 
has been called a topological form of nematic order by 
Lammert et al. 20] in a generalized Z2 gauge theory of 
nematic order. Also there rotational (Lorentz-) invari- 
ance is restored even though there is no condensate of 
disclinations, and the Burgers vectors are disordered by 
fluctuations (see also ref. jlal)- The dual of this phase is 
Coulomb-like. 

The above description of defects was formulated in 
what has been named tangential approximation to the 
Euclidean group 0, in which the discrete rotations are 
treated as if they took place in the tangential place with 
arbitrary real Franck vectors SI' [see Eq. (|1.4|l ]. In a 
more accurate formulation, the nonabelian nature of the 
rotations and the quantization of Sl^ must be taken into 
account. Their discreteness is certainly remembered in 
the nematic phase, even if the directions of the Burg- 
ers vectors become diordered (see also the discussion 
in Ref. ^^). This implies that there are elements of 
quantized curvature fluctuating in spacetime. Fortu- 
nately, this does not introcduce any unphasical results 
at presently accessible length scales since these fluctua- 
tions mainly renormalize the basic curvature energy in 
(|1.20l) . as discussed before. 

The generalization to four Euclidean spacetime dimen- 
sions changes mainly the geometry of the defects. In four 
dimensions, they become world sheets, and a second- 
quantized disorder field description of surfaces has not 
yet been found. But the approximation of representing a 
sum over dislocation surfaces in the proliferated phase as 
an integral as in Eq. (|1.18|) will remain valid, so that the 
above line of arguments will survive, this being a natu- 
ral generalization of the Meissner-Higgs mechanism. The 
disclination sources of curvature will be world sheets, as 
an attractive feature for string theorists. However, the 
high-energy properties will be completely different. On 
the one hand, these surfaces behave nonrelativistically 
as the energies approach the Planck scale, on the other 
hand they will not have the characteristic multi-Planck 
recurrences of the common strings. Although the lat- 
ter property may never be verified in the laboratory, the 
deviations from relativity at high energies or short dis- 
tances may come into experimentalists reach in the not 
too distant future. 

Note that our model has automatically a vanishing cos- 
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mological constant. Since the atoms in the crystal are in ilar to that given by Volovik '2'| with his hehum droplet 
equilibrium, the pressure is zero. This explanation is sim- analogies. 
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